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, Abstract. For confining potentials of the form q{r) = r + p(r), where p(r) decays rapidly and is smooth for 

■ r > 0, it is proved that q(r) can be uniquely recovered from the data , Sj}vj=i,2,3.... • Here Ej are energies 

' of bound states and Sj are the values n^ (0), where Uj{r) are the normalized eigenfunctions, u^dr = 0. An 

algorithm is given for finding g(r) from the knowledge of few first data, corresponding to 1 < j < J assuming 
^jQ, that the rest of the data are the same as for qo{r) := r. 

< 

1. Introduction. 

The problem discussed in this paper is: to what extent does the spectrum of a quarkonium system 
together with other experimental data determines the interquark potential? This problem was discussed 
in [1], where one can find further references. The method given in [1] for solving this problem is this: 
, one has few scattering data Ej, Sj, which will be defined precisely later, one constructs using the known 

CD ■ results of inverse scattering theory a Bargmann potential with the same scattering data and considers 

this a solution to the problem. This approach is wrong because the scattering theory is applicable to 
the potentials which tend to zero at infinity, while our confining potentials grow to infinity at infinity 
I and no Bargmann potential can approximate a confining potential on the whole semiaxis (0, cxd). The 

aim of this paper is to give an algorithm which is consistent and yields a solution to the above problem. 
The algorithm is based on the well-known Gelfand-Levitan procedure [2]- [4]. 
Let us formulate the problem precisely. Consider the Schroedinger equation 

-V^jj -I- q{r)i:j = EjiPj in R^ (1.1) 
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^ I where q{r) is a real-valued spherically symmetric potential, r :— \x\,x £ R^, 

^ , q(r) = r +p(r), p(r) = o(l) as r — s- oo. (1-2) 

•1—1 , 

^ . The functions ipji^): \\^j\\L^{R^) = 1; are the bound states, Ej are the energies of these states. We define 

$— ( ' ''^j(^) •= ''V'j('')i which corresponds to s- waves, and consider the resulting equation for uj: 

Luj := -u'- +q{r)uj = EjUj, r > 0,Uj{0) = 0, ||uj||i2(o,oo) = 1- (1-3) 

One can measure the energies Ej of the bound states and the quantities Sj = u'^ (0) experimentally. 
Therefore the following inverse problem (IP) is of interest: 
(IP): given: 

{-B„s,}v,=i,2,... (1.4) 

can one recover p{r)7 

In [1] this question was considered but the approach in [1] is inconsistent and no exact results are 
obtained. The inconsistency of the approach in [1] is the following: on the one hand [1] uses the inverse 
scattering theory which is applicable only to the potentials decaying sufficiently rapidly at infinity, on 
the other hand, [1] is concerned with potentials which grow to infinity as r —^ +00. It is nevertherless 
of some interest that numerical results in [1] seem to give some approximation of the potentials in a 
neighborhood of the origin. 

Here we present a rigorous approach to the problem considered in [1] and prove the following result: 
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Theorem 1. IP has at most one solution and the potential q{r) can be reconstructed from data (1-4) 
algorithmically. 

The reconstruction algorithm is based on the weU known Gelfand-Levitan procedure for the recon- 
struction of q{x) from the spectral function. We show that the data (1.4) allow one to write the spectral 
function of the selfadjoint in L^(0, oo) operator L defined by the differential expression (1.3) and the 
boundary condition (1.3) at zero. 

In section 2 proofs are given and the recovery procedure is described. 

Since in experiments one has only finitely many data {Ej, Sj}i<j<j, the question arises: 

how does one use these data for the recovery of the potential? 

We give the following recipe: the unknown confining potential is assumed to be of the form (1.2) and 
it is assumed that for j > J the data {Ej, Sj}j>j for this potential are the same as for the unperturbed 
potential qo{r) = r. In this case an easy algorithm is given for finding q{r). 

This algorithm is described in section 3. 

II. Proofs. 

We prove Theorem 1 by reducing (IP) to the well-studied and solved problem of recovery of q{r) from 
the spectral function [2], [3]. 

Let us recall that the selfadjoint operator L has discrete spectrum since q(r) +oo. The formula for 
the number of eigenvalues (energies of the bound states), not exceeding A, is known: 



Vl:=A^(A)~-/ [X-q{rt^dr. 



Ej<\ ^ Jq(r)<\ 

This formula yields, under the assumption g(r) ~ r as r — > oo, the following asymptotics of the eigen- 
values: 

^ as J ^+00. 

The spectral function p(A) of the operator L is defined by the formula 



(2.1) 



where aj are the normalizing constants: 



rj{r)dr. (2.2) 

Here (/>j(r) := (j){r, Ej) and (/)(r, E) is the unique solution of the problem: 

L(t) := -(j}" + qir)<f) = Ecp, r > 0, (?i(0, E) = 0, (?i'(0, E) = 1. (2.3) 

If E = Ej, then (pj = (j)(r,Ej) e L^(0,oo). The function <p{r,E) is the unique solution to the Volterra 
integral equation: 

,, sm{\/Er) r sin] VE(r — y)] ,„ 

<P{r,E)= ^r-' + ^ )- ^'^ q{y)ct>{y,E)dy. (2.4) 

vE Jo VE 

For any fixed r the function 4> is an entire function of E of order |, that is, < cexp(c|-E'|^/^), 
where c denotes various positive constants. At E = Ej, where Ej are the eigenvalues of (1.3), one has 
(f>{r, Ej) := 4>j E L^(0, oo). In fact, if q{r) ~ cr°, o > 0, then \(j)j\ < cexp(— 7r) for some 7 > 0. 
Let us relate aj and Sj. From (2.3) with E = Ej and from (1.3), it follows that 
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Therefore ^ 

:= ll'?^jlli,2(0,oo) = 3- (2-6) 

Thus data (1.4) define uniquely the spectral function of the operator L by the formula: 

p(A) := ^ si (2.7) 

Ej<X 

Given p(A), one can use the Gelfand-Levitan (GL) method for recovery of q{r) [2], [3]. According to this 
method, define 

a{X) := p(A) - po(A), (2.8) 

where /9o(A) is the spectral function of the unperturbed problem, which in our case is the problem with 
q{r) = r, then set 

poo 

L{x,y):= Mx,^)My,^)dcr{X), (2.9) 



where (j)o{x,X) arc the cigcnfimctions of the unperturbed problem (2.3) with q{r) = r, and solve the 
second kind Fredholm integral equation for the kernel K{x, y): 



K{x, y)+ f K{x, t)L{t, y)dt = -L{x, y), < y < x. (2.10) 
Jo 

The kernel L{x,y) in equation (2.10) is given by formula (2.9). If K{x,y) solves (2.10), then 

p^r) = 2^^p^, r > 0. (2.11) 
dr 

3. An algorithm for recovery of a confining potential from few experimental data. 

Let us describe the algorithm we propose for recovery of the function q{x) from few experimental data 
{Ej, Sj}i<j<,7. Denote by {E^j ■ •Sj}i<i<,/ the data corresponding to qo := r. These data are known and 
the corresponding eigenfunctions (1.3) can be expressed in terms of Airy function Ai(r), which solves the 
equation w" — rw = and decays at +oo, see [5]. The spectral function of the operator Lq corresponding 
to g = g'o := ^ is 

po(A) := ^ (4)^. (3.1) 

Define 

p(A):=po(A)+a(A), (3.2) 

■■= E 4 - E (^°)'' (3-3) 

Ej<X E°<\ 

and 

J J 
L{x,y) ■.= J2s'^^{x,Ej)cj>{y,Ej) - J2{s'ifM^)My), (3-4) 

(p{x,E) can be obtained by solving the Volterra equation (2.5) with q{r) = qQ{r) := r and represented 
in the form: 

cPix, E) = ' + K{x, y)^^^dy, (3.5) 

where K{x,y) is the transformation kernel correponding to the potential q{r) = qo(r) := r, and (pj are 
the eigenfunctions of the unperturbed problem: 

+ r^j = Ej(t>j r > 0, (t>j{Q) = 0, <^;(0) = 1. (3.6) 
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Note that for E ^ E'^ the functions (3.5) do not belong to L'^{0,oo), but ((){0,E) = 0. Wc denoted 
in this section the eigenfunctions of the unperturbed problem by <j)j rather than cpoj for simphcity of 
notations, since the eigenfunctions of the perturbed problem are not used in this section. One has: 
(f>j{r) = CjAi{r — Ej), where cj = [Ai' {— Ej)]~^ , Ej>Ois the j— th positive root if the equation 
Ai{—E) = and, by formula (2.6), one has = [cj Ai'^{r — Ej)dr]^^/^ . These formulas make the 
calculation of 4>j{x), Ej and s° easy since the tables of Airy functions are available [5]. 
The equation analogous to (2.10) is: 

2J px 27 

K{x,y) + J2^j'^jiy) K{x,t)^j{t)dt = -J2cj^ji^)'^jiy)^ (3-7) 

where *j(t) := (/)(t, Ej), cj = s|, 1 < j < J, and 5'j(t) = (/)j(t), c-j = (s°)2, J+ 1 <j< 2J. Equation (3.7) 
has degenerate kernel and therefore can be reduced to a linear algebraic system. 
If K{x,y) is found from (3.7), then 

p(r) = 2-^K{r,r), q{r) = r+p{r). (3.8) 

Equation (2.10) and, in particular (3.7), is uniquely solvable by the Fredholm alternative: the homoge- 
neous version of (2.10) has only the trivial solution. Indeed, if h + L{t,y)h{t)dt = 0,0 < y < x, 
then + J^^\h\'^[dp{X) - po{\)] = 0, so that, by Parseval equality, J^^\h\^dpi\) = 0. Here 

h := h{t)(f>{t, X)dt, where (j){t. A) are defined by (3.5). This implies that h{Ej) = for all j = 1, 2, .... 
Since h{X) is an entire function of exponential type < x, and since the density of the sequence Ej is 
infinite, because Ej = 0{j^^^), as was shown in the beginning of section II, it follows that h = and 
consequently h{t) = 0, as claimed. 

In conclusion consider the case when Ej = Ej, Sj = s° for all j > I, and {Eq, sq} is the new eigenvalue, 
£^0 < -£'1) with the corresponding data sq. In this case L{t,y) = SQ(f)o{t, Eo)(f)o{y, Eq) , so that equation 
(2.10) takes the form 

K{x,y) + slMy) I K{x,t)Mt,Eo)dt = -sl(j)o{x,Eo)My^Eo). 
Jo 

Thus, one gets: 



drl + sl£<j)l{t,Eo)dt' 
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